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Abstract 

We consider the Standard Model with an arbitrary number nn of Higgs 
doublets and enlarge the lepton sector by adding to each lepton family £ a 
right-handed neutrino singlet uir. We assume that all the Yukawa- coupling 
matrices are diagonal, but the Majorana mass matrix Mji of the right-handed 
neutrino singlets is an arbitrary symmetric matrix, thereby introducing an 
explicit but soft violation of all lepton numbers. We investigate lepton- flavor- 
violating processes within this model. We pay particular attention to the 
large-m/j behavior of the amplitudes for these processes, where ruR is the 
order of magnitude of the matrix elements of Mr. While the amplitudes for 
processes like M 7 and Z drop as ^/m?R for arbitrary uh, 

processes like ji^e^e" and ^~ — > e~e+e~ obey this power law only for 

uh = 1. For uh > 2, on the contrary, those amplitudes do not fall off when 
rriR increases, rather they converge towards constants. This "non-decoupling" 
of the right-handed scale occurs because of the sub-process £~ — > £'~ S^* , 
where is a neutral scalar which subsequently decays to e~^e~. That sub- 
process has a contribution from charged-scalar exchange which, for uh > 2, 
does not decrease when rriR tends to infinity. We also perform a general study 
of the non-decoupling and argue that, at the one-loop level, after performing 
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the limit rriR — > oo and after removing the i/ji from the Lagrangian, our model 
becomes a normal multi-Higgs-doublet Standard Model with loop-suppressed 
flavor-changing Yukawa couplings. Finally, we show that in our model the 
branching ratios of all lepton-fiavor-changing processes are several orders of 
magnitude smaller than present experimental limits, if one makes the usual 
assumptions about the mass scales in the seesaw mechanism. 

PACS numbers: 14.60.St, 13.35.-r, 11.30.-j 
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I. INTRODUCTION 



Recent experimental evidence strongly suggests that neutrinos mix and are, therefore, 
massive • This raises the question of why are the neutrino masses so much smaller than 
the masses of charged fermions. A simple answer to this question is provided by the seesaw 
mechanism In a seesaw model there are two mass scales: 

m^), the scale of the Dirac mass terms linking the (known) left-handed neutrinos to (new) 
right-handed neutrinos; 

mn, the scale of the Majorana mass terms of the right-handed neutrinos. 

When mji ^ mr) the mixing of the two types of neutrinos gets suppressed by mjy/mR and 
the left-handed neutrinos acquire Majorana masses of order m'jj/mR, hence much smaller 
than the Dirac mass terms. Now, m/j ^ m/j is natural, since the Majorana mass terms of 
the right-handed neutrinos are gauge-invariant; therefore, they do not need to be of order 
of the Fermi scale mp of the spontaneous breaking of the electroweak symmetry. 

In the context of the seesaw mechanism an interesting option consists in having lepton- 
flavor symmetries which are respected by the Dirac mass terms, but broken by the Majorana 
mass terms of the right-handed neutrinos 0,0] • This option arises because the Dirac mass 
terms originate in Yukawa couplings of the leptons to scalar doublets, which are hard (di- 
mension four), while the Majorana mass terms are soft (dimension three). Lepton- flavor 
breaking thus becomes soft and, in fact, an unsuppressed reflection at low energies of some 
physics at ultra-high energies. This attractive hypothesis has the advantage that it allows 
one to construct simple models which explain the apparent maximal mixing of atmospheric 
neutrinos, solar neutrinos, or both simultaneously J^. 

Once it is accepted that lepton-fiavor breaking is soft, there is nothing against introducing 
many scalar doublets with Yukawa couplings to the leptons, since lepton-fiavor-changing 
neutral interactions are automatically absent at tree level from those couplings. The question 
then arises of knowing whether lepton-fiavor-violating decays, which arise at loop level, are 
suppressed by some powers of m£,/mR or mp/mR, or not. Moreover, one would like to 
identify the effective field theory at low scale which corresponds to the limit ttir — oo. 

In this paper we try and answer the questions above by computing the lepton-fiavor- 
violating decay r~ —>■ ^~e^e~ in the context of a seesaw model with an arbitrary number 
of scalar doublets and with softly-broken lepton numbers. We take this tau decay as a 
concrete example for the study of general features of our model. Since we compute the full 
one-loop decay amplitude for r~ — > ;U~e"'"e~, we simultaneously also have the amplitudes 
for r~ At~7, Z t"*"/!", and — > r+/x~ at our disposal. (5"° denotes a physical neutral 
scalar.) Evidently, the changes from tau decays to muon decays can be performed trivially in 
our results. There are in the literature a number of analogous computations (see Refs. [§-[1^ 
and citations therein), yet our work is different for the following reasons: 

i. In previous works the seesaw mechanism is considered in the context of the Standard 
Model, i.e. with only one Higgs doublet. There are then no physical charged scalars, 
and the neutral scalar — the Higgs particle — has been neglected because its Yukawa 
couplings are suppressed by the smallness of the charged-lepton masses. In the present 
work we consider charged scalars in the loops, and also neutral scalars 5*^ in the process 
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r- ^ fi-Sr /i"e+e- . (1.1) 

ii. In previous works all external momenta have been set to zero. In this paper we 
give exact expressions for non-zero external momenta. Since the masses of the light 
neutrinos are much smaller than the masses of the charged leptons, it does not seem 
justified to treat the former exactly while neglecting the latter. 

iii. We study how the various contributions to the decay amplitude behave as functions 
of m^. We demonstrate that the contributions previously computed are proportional 
to l/mj^, and thus negligible for sufficiently high mn, while on the other hand some 
contributions to the process ( |1.1| ) remain unsuppressed. 

In order to perform our computation in the context of a general multi-scalar-doublet 
model, we had to extend the formalism previously developed for the scalar particles in that 
model This formalism is presented in detail in Appendix it may be useful for other 



computations in that general model. We also took a close look at the one-loop renormaliza- 
tion of flavor-changing interactions; we show in Appendix ^ that the fermion wave-function 
renormalization, including fermion mixing, does not introduce any contributions to flavor- 
changing decays beyond those given by diagrams with flavor- changing self-energies in the 
external fermion legs [|15| . 



This paper is organized as follows. In Section |T| we discuss the leptonic couplings and the 
necessary formulas concerning the seesaw mechanism. Section |ITT| treats some notation for 



the process r~ — > /i~e"'"e~ , our example decay. Section |^ deals with the orders of magnitude 
in our model and Section |V| introduces conventions and sub-processes for the example decay. 
Sections ^I|, |VII| , |VIII| and |X| describe, respectively, the photon, Z, neutral-scalar, and 



box-diagram sub-processes of r~ fi~e^e~ . In Section ^ we discuss the limit of infinite 
right-handed scale. Section ^ presents decay rates for the example decay and for other 
fiavor- changing decays whose amplitudes have been implicitly calculated in Sections [V]| - 
VIII| . The conclusions are found in Section pCIl . 



II. THE LEPTONIC COUPLINGS 

A. General seesaw framework 

We consider an extension of the standard model with three families and three right- 
handed neutrinos. We label the latter with family lepton numbers: UeR, i^^r-, and Vt-^. At 
the moment we do not assume conservation of the family lepton numbers, thus this labelling 
has no physical content. The Yukawa Lagrangian of the leptons is 



+ H.C. 



k=l i,£'=e,fj.,T \ / 

(2.1) 

The mass matrix of the charged leptons and the Dirac neutrino mass matrix Md are 
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Me^^y^vlTk and Mo = i=V^;,Afc, (2.2) 

respectively. Without loss of generality, we assume iW^ to be diagonal with real and positive 
diagonal elements: = diag (mg, m^, m,-). The mass terms for the neutrinos are 

- VrMdi^l - ^ vrCMrVI + H.c. , (2.3) 

where C is the charge-conjugation matrix and Mr is non-singular and symmetric. 

The left- and right-handed neutrinos are written as linear superpositions of six physical 
Major ana neutrino fields xi'- 

^iL = ^{UL)^i-iLXi and l^^R = ^{UR)^^-tRXi, (2.4) 

i i 

where 7^ = (1 — 75) /2 and 7r = (1 -|- 75) /2 are the projectors of chirality. The fields Xi 
satisfy Xi = Cxf- The matrix 



TT* 



is 6 X 6 unitary; therefore, 



and 



C/lC/1 = 13x3, (2.6) 

Wi = l3x3, (2.7) 

UlU^^ 03^3, (2.8) 



UlUL + UlU*j,^le^Q. (2.9) 



U is defined in such a way that 



^^iH ) W = m = diag(mi,m2,...,m6), (2.10) 

with real and non- negative m^. Therefore, 

f/2mf/l = 0,3x3, (2.11) 

URmUl = MD, (2.12) 

URmUl = MR. (2.13) 

The charged-current Lagrangian is 

= -4 E K" ^^^^^^ ^"^'^^^^ + (^^) ^'^'^^^] • 
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The interaction of the Z boson with the leptons is given by 



Xj 



(2.15) 
(2.16) 



When extracting the vertex from Eq. (|2.15|) one must multiply by a factor 2, since the 
neutrinos are Majorana fields. 

The Yukawa couphngs of the charged scalars S*^ to the leptons are written in the following 
general notation: 

= {San{Ra)ulR - {La),M X^ + S+X^ [K),,7L " (4),,7r] ^} • (2-17) 

The notation for the scalar sector, and the precise meaning of the nj^- vectors a, are explained 
in Appendix |A 1|. One has 



Ra = AIUr , withA„ = J2 (^k^k 

k 

La = TaUL , withFa = ^ a^Ffc . 



(2.18) 
(2.19) 



Similarly, the Yukawa couplings of the neutral scalars (see again Appendix [A 1| ) to the 
leptons are written 



^ b L i,j '- 

+ E4(r,),,,7. + (rt) 7« 



(2.20) 

(2.21) 



with Th = "^k^k^k- The matrices Fb are symmetric: 

with Afc = ^ bkAk . 



(2.22) 



When extracting the Feynman rule for the vertex with the neutrinos from Eq. ( ^.2U| ), one 
should multiply by a factor 2, since the Xi are self-conjugate fields. 

In the case of the charged Goldstone bosons = S^^, one has (see Appendix A. 3) 
iaw)k = (2mvi/)- Therefore, 



A, 
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aw 



r 
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V2mw 



Md 



(2.23) 

(2.24) 
(2.25) 
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and 



R. 
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■aw 



-'aw 



V2mw 
9 



Uirh, 



(2.26) 
(2.27) 



B. Our model 



We assume that the Yukawa- couphng matrices Tk and are simultaneously diagonal 
0. Therefore, the matrices Fq, A^, Ffe, and A;, are all diagonal. The neutrino Dirac mass 
matrix M^j is also diagonal, hence 

ULfhU^R = Ml and Ultti^uI = M^Md are diagonal. (2.28) 

Now the labelling of the neutrino fields according to family lepton numbers acquires a well- 
defined meaning. 

Diagonal Yukawa-coupling matrices are achieved by assuming invariance of the Yukawa 
Lagrangian under U (1) (a = e, /i, r), the groups associated with conservation of the lepton 
number for each lepton family. Since the gauge part of the Standard-Model Lagrangian 
is invariant under these U{1) symmetries anyway, and since the scalar doublets do not 
transform under these U{1) groups, the only place where these lepton-number symmetries 
are violated is the Majorana mass term of the right-handed singlets in Eq. ( |2.3| ). Since the 
mass term is an operator of dimension three, this violation is soft [0. Hence, the one- loop 
amplitudes for lepton-flavor-violating processes must be finite. 



III. NOTATION FOR THE PROCESS 

The flavor-changing decay that we want to study is 

(Pi) (P2) e+ (ps) e~ (pa) . (3.1) 



Clearly, 



We denote 





pI 


= ml 




pI 


= ^1 


pI 


= pI 


= ml 



(3.2) 
(3.3) 
(3.4) 



q» = p» -p» = p» +p2 , 
r"" = Pi — p2 = P2 + Ps 5 



and 
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Pi = • (3.6) 



Then, 

Ami < (f < i^r - m^f , 



{m^ + rUe) < r < (m-r- — rrie) . 

The amphtude M for the process of Eq. involves the Dirac spinors = {P2), 

Ur = Ur (pi), Ue = Ue (^4), and Ve = fe (ps)- These spinors satisfy 

«/.^2 = rrif^u^ , (3.8) 
]^iUr = rrirUr , (3.9) 

Ue^4. = meUe, (3.10) 
^3Ve = -meVe. (3.11) 

In our calculations we shall often need the the following product of elements of Ul: 

^^ = (Ul)^, [Ul)^^ , satisfying ^x, = , (3.12) 

i 

where the second relation is a consequence of unitarity — see Eq. ( |2.6| ). Furthermore, we shall 
make use of dimensional regularization, evaluating integrals in a space-time of dimension d. 
Therefore, we define 

d = A-2e and Vk = d'^k/ {2nf , (3.13) 

where the latter expression is an abbreviation for integration over the momentum k. Even- 
tually, we shall take the limit e ^ 0. Then, in some integrals the divergent constant 

/C = e-^ - 7 + ln(47r) (3.14) 

will appear (7 is Euler's constant). The reason why /C always drops out when calculating 
the amplitude for the process ( p.l| ) in our model will be discussed in detail. 



IV. ORDERS OF MAGNITUDE 

We assume that the matrix elements of Md are of order mo and the square roots of the 
eigenvalues of Mj^Mpt are of order mpt, with m^) <^ m^j. Then, 

mi ~ and {Ur)„. ~ for i 

mn mn 

mi ~ mR and {ULjoj ~ — ^ for i 

There is also the order of magnitude, which we may call m^, of the charged- lepton masses. 
This may be taken as either m^- ~ 2 GeV, or m^ ~ 10~^ GeV, or mg ~ 5 x 10~^ GeV. 
Because of Eqs. (p. 7]), q"^ m'j. For simplicity we shall identify with mn- 



1,2,3; 



4,5,6. 



(4.1) 
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Finally, there is the Fermi scale mp, in between 10 GeV and 1 TeV. The W mass mw, 
the charged-scalar masses nia, and the neutral-scalar masses are all taken to be of order 
mp. The scalar-potential couplings Caa'b in Eq- ( |A30| ) are also assumed to be of order mp. 



The overall hierarchy of mass scales is thus mp, -C mp <^ m^. 

When we have information about the masses of the light neutrinos we can estimate the 
order of magnitude of mR via the seesaw relation ||^ in the first line of Eq. (|4.1| ). If we take 
the light-neutrino masses to be of the order of y^Am^m ~ 0.05 eV, where Am\^^ is the 
neutrino mass-squared difference relevant for atmospheric-neutrino oscillations , and if we 
assume that mp) ~ or m,-, then we obtain mp^ ~ 10^ -h 10^^ GeV. Thus we might regard 
mpi ~ 10^'^ GeV as a typical order of magnitude of the right-handed scale, keeping in mind 
however that mp might deviate several orders of magnitude from this value. 

It is important to notice that, with the convention m^ = mp,, the matrices Ra^^ and 
Layy in Eqs. ( p.26|) and (|2.27| ), giving the Yukawa couplings of the charged Goldstone boson. 



are of the same orders of magnitude as the corresponding general matrices and La in 
Eqs. (|2.19| ) and (p.l8|) , provided we make the assumption 

Tj ,Aj mp/mp; (4.2) 

this is a natural relation in view of = mp. This means that the factors m|,/m|. suppress- 
ing some charged- Goldstone-boson contributions are exactly the same as those suppressing 
the corresponding, and more general, charged-scalar contributions. 

In the result for the momentum integrals large logarithms, like for instance \n{mp/mp), 
arise. We shall not take into account such logarithms in our estimates of orders of magnitude. 



V. CONVENTIONS AND SUB-PROCESSES 

We shall compute the process ( ^.11 ) using the conventions and vertices given in Ref. [IB . 
In the Z^-G^ sector we use the unitary gauge, thus discarding G°. On the contrary, in 
the W^-G^ sector we use Feynman gauge. This means that the propagator of is 
—igpa/{k'^ — m^). The charged- Goldstone-boson contributions will usually be taken into 
account together with the general charged-scalar contributions; the sums over the charged 
scalars 5*^^ will not exclude the charged Goldstone boson G^ = . We remind that, in 
Feynman gauge, = mw 

The process may proceed via box diagrams or through one of the three following 
sub-processes. In the sub-process with amplitude M4 the initial lepton r~ decays into the 
final lepton together with a virtual photon A with momentum g; the photon later decays 
into the e^e~ pair in the final state. The corresponding amplitude may be written 

Ma = ^eUe^ipVe, (5.1) 

where e is the positron charge and is the amplitude for r~ (pi) /i~ (^2) 7 (q'), the r~ 
and the /i~ being on mass shell while the photon is off mass shell. Current conservation, 
i.e. Eqs. ( |3.10| ) and ( p.ll| ), implies that one may discard all terms proportional to from 
M^. The sub-process with amplitude Mz is analogous to the previous one with the virtual 
photon substituted by a virtual Z boson. Its amplitude may be written 
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Finally, there is the sub-process with amplitude Mh, in which the r decays into /i together 
with an off-shell neutral scalar with momentum q. The corresponding amplitude is 



Mh = —= Me 



ee 



(r,L7L+ rn v,- (5.3) 



One must sum over all physical neutral scalars h ^hz- 

Besides these three sub-processes, there are also box diagrams, which are all finite, to be 
considered in Section IIX. 



No one-loop diagram with a neutral scalar in the loop can contribute to the process 

3) 

The infinities in the amplitude for the process ( p.l|) cancel for the following reasons: 

A. Conservation of the electromagnetic current; 

B. unitarity of the diagonalization matrix W; 

C. flavor- diagonal Yukawa-coupling matrices. 

Item A is independent both of our model and of the seesaw mechanism, and applies to 
the photon sub-process. Concerning item B, the relations ( |2.6|) and (|2.8| ) are relevant, the 
first one in the form of Eq. ( p.l2| ). Only item C is directly connected with our model and, 
clearly, it plays a role only in charged- scalar exchange. More generally, items B and C are 
responsible for the cancellation of all terms independent of the neutrino masses mj. 



VI. r 



A. Computation 

Let us work out in detail. First consider the transition T~{p) —>■ fi^{p), effected 
either by exchange or by charged-scalar exchange. Call the corresponding amplitude 
— zS(p). Then, there are the following two contributions to M^: 



We obtain for T,{p) 
- i^{p) = 

where 



ml — m?, 



-eu. 



-iT. {p2)\ — — YUr . 



— 



vkY — 



2 {k-pf-rri^ 



E 



(k-p)- ml 



(6.1) 
(6.2) 



(6.3) 
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-m, (L,)^, {Ri),,lL - (Ra)^^ (4),.7i? • (6-4) 



The first term in the right-hand-side of Eq. ( |6.3|) is the l^^-exchange contribution. The 
second term automatically includes, in the sum over a, the contribution of the charged 
Goldstone boson — which is obtained by setting rriayy = mw and by using Eqs. (|2.26| ) and 
( p.27| ) for the matrices Ra^^r and Law- 
There is one more contribution to MjJ, namely from diagrams in which the photon 
attaches either to the or to the charged scalar in the loop. There is no W^S'^'y vertex 
except in the case 5 J = G^; also — see Appendix |A4| — the vertex W^G'^'j has a factor 
mw which cancels the denominator in the Yukawa coupling of . One obtains the 
amplitude 



/^^ E ^ ( n n - 2) - P+r^lL - 4fc ■ p+YIl + ^f+¥lL 

{Y^-m IR + m, Q^Y - Yfi) 1l] + Yl I ■ (6.5) 

a J 

We have used the shorthands 

D,w = {k- p,f - and Aa = {k - p,f - ml, for l = 1,2. (6.6) 

The sum over a in Eq. ( |6.5| ) includes the contribution of the vertex G^G~'~f. The other terms 
in the right-hand-side of that equation are the contributions from the vertex W~^W~'j and 
from the vertices W^G'^j. 

Thus, = M^, + + M^3. 

One writes the Feynman integrals symbolically: 

'^kj^ = a\, [Vk-^=b\pl; (6.7) 

T^k \ = 4 . / J][ = bipl ; (6.8) 
DiD2w J DiD2w 



Vk ^ = , fvk = d^p\ + , 

DiDiwD2w J DiDiwD2w 



(6.9) 



In a\,...,u\ the index i indicates the dependence on mf. The coefficients u*, a\, a\, b\, 
and b\ diverge when (i — ^ 4, while all other coefficients are finite. However, the following 
relations hold: 
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a - a; 



<f (4 - 4) + {^1 - ml) , 



(6.10) 

b\ - bl = {d\ - 4) + {ml - ml) {k{ + k\) , (6.11) 
mlU, - mlbi = {mld\ - m^^) + [ml - ml) {2u' - q'f + mlk\ + mlk\) , 

(6.12) 



where 



k\ = d\ -\- — c\ , 
k2 — ^2 -\- f ~ C2 , 

/C3 = C]^ + C2 — 0-3 . 



(6.13) 
(6.14) 
(6.15) 



One may write relations similar to Eqs. (|6.7| )- (|6.15| ) for the integrals which have D^a in- 
stead of D^w (for = 1)2), then we use the notation Oi", . . . in order to indicate the 
dependences on m^ and on m^. 

Using Eqs. ( |6.10| )- p.l2| ) one finds that is finite and that, moreover, it respects gauge 
invariance, i.e. qpM'^ = 0: 

M^ = e^Up {iaf'^qx {au 7l + oiRi 1r) 



ml — m?, 



q'^lR 2 



q^ 



mi — mt 



■ Y {t^iiIr + "^t7l) 



(6.16) 



One obtains (for d = 4) 

OiLi = g'^Xim^ {kl + c\ + 4) 



ttffi = g^Ximr {k\ + Ci + 4) 

+ ^ [m. (i?J^, {Ri)^^ k^ + (LJ^, (lI)^^ - (i?J^, (lI) A^^' 

a 

/3Li = g'^x-im^ [kl + ^ + 4) 

+ 5^ [m, (i?J^^ [Ri)^^ kr - (LJ^, (Lt)^^ fc^ + (^a),. A4"^' 



(6.17) 



(6.18) 



= -g'^ximr {kl + Ci + 4) 

+ 5^ [-m. {R^)^^ (Rl)^^ kr + m, (L„)^^ (4)^^ + (i?J^^ (4)^^ Ac^-'^' 



(6.19) 



(6.20) 
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where 



4 = r-4, (6.21) 

4 = f-4, (6.22) 



and similarly for the coefficients ^4'^ and fcg". Also, Ac*'"^ = c<^ — . 



B. Order of magnitude 

Integrating over the momentum k one finds, from Eqs. ( |6.9|) , 

>1 /-l-x 



(03,4,4,4,4,/*) = j dx j dy—^{l,y,x,y^,x'^,xy), (6.23) 

where 

Af^ = {x + y) + {1 — X — y) {ml — xm^^ — ym^) — xyg^ . (6.24) 

Let us denote by a generic coefficient 4, 425 '^125 similarly (7*'" for 4''^; 

etc. Clearly, is a function of mf, with gf* ~ rrip^ for i = 1,2,3 and g^ ~ for 
i = 4,5, 6. Moreover, for i = 1,2, 3, is independent of up to corrections of order 
mf/rn^ ~ '^i)/ (^r^f)- The same holds for the (7*'", of course. 

Considering Eqs. (|6.16| )- (|6.2CI| ) in detail, we see that involves terms of five types: 

1. g^mjY^.Xig'; 

2. rnjEARa),,{Ri),,g''''; 

4. meE.mdRa),,{Li)^^g^'^; 

5. m,E.'^.(^a),,(i?i),,^7*'^ 

We remind the reader that we set m£,. Also, from Eqs. ( |2.18D and ( |2.19| ), we see that 

La is Ul and Ra is Ur except for diagonal Yukawa- coupling matrices. 

Terms of types 1 and 3 are similar. They are proportional to {UL)^i{U'[)ir, which is of 
order 1 for i = 1,2,3 and of order m'jj/mj^ for i = 4,5,6. As we have seen, g^ ~ rrip^ is 
almost independent of i for i = 1,2, 3. One may therefore use 

3 62 
m 



Esr^ "I'D 



(6-25) 



mo 

i=l i=i ^ 

to estimate mj X^Li ^i9^ t)e of order m^/ (m^m|,). For i = 4, 5, 6 one has g^ ~ m^^ and 
Xj ~ m'j^/m'j^, hence terms of type 1 or 3 are of order mj^/m'^. 

Terms of type 2 are proportional to {Uji)^i{Ul^iT. Therefore, for i = 1,2,3 they are of 
order m^,/ [m\my) and for z = 4, 5, 6 they are of order m\)/m\. 
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Finally, terms of types 4 and 5 always include a suppression mo/mji from the mixing 
matrices Ul and Ur. Taking into account also m^, m^, and g^'"", those terms are of order 
m\)/ {m\rn?p) for i = 1, 2, 3, or mjj/m'j^ for i = 4,5, 6. 

In summary, we find that the amplitude for the vertex fi^A* is suppressed by 

m^/m^. With mo ~ 1 GeV and mji ~ 10^'^ GeV one has a suppression factor 10^^° in 
the amplitude of r ^ fi'j. Note that terms of type 2-5 contain a product of two Yukawa 
coupling constants. If we are more specific and assume the relation ( [4 .21 ) for the Yukawa 
couplings, then there is a suppression factor of order m|)/(m|,m|n) for all of them and, in 
addition, a (IGvr^)"^ from the loop integration. Terms of type 1, on the other hand, have an 
extra factor g"^. 



VII. T 



A. Graphs in which the Z attaches to charged particles 

The vertex fi~Z* has three contributions M^^, M^2, and M^g analogous to the 

corresponding contributions to fi^A*. Comparing them one easily concludes that 



'Z2 



(7.1) 



with 



— mi 



+u. 



g^Xi 



[(2 -d){k- p+fJ^^L + 4A; ■ p+YlL - ^pif^lL 



DiDmrD2w 

+mr}iYlR + m^Y¥lL - 2m^^YlL] Mr • (7.2) 
We compute by using the same method as in Eqs. (|6.7| )- (|6.22|) , arriving at the result 

T' = - {u,lRUr)Y,9'x,mr [2p^ [kl + c\) - [k^ + c\)] 

i 

- (^M7L^^r) 5^ (7'a;,m^ [2p^ (4 + 4) + q" {k^ + 4)] 

i 

^ {Lai + Ran) Ur , (7.3) 



where 



L, = g^Xi [{2 -d)u'- 2m24 + [2ml + ^ - 4 + [ml + 2mJ - q^) 4] 



R 



i = g Xim^m^ 



C-^ ~\~ k-^ H~ 



mt — mf 



(4 - d\) 



-kY 



kr + 



mt — m 



[di - d{ ) 



(7.4) 
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o i.a I 2 ci.a 2 ; i,a, 2 ; «,a i ^ / 2 2 j«,ct\ 

-2m • + g / • - m^k{ - m + — [m - m^d{ ) 



+ (La),, {Li) 



ml — ml 



( i,a i,a\ 
(Cl - C2' j 



(7.5) 



Let us consider the order of magnitude of T^. We realize at a glance that most of its 
terms are of one of the five types already considered in Subsection |V1 B| . In Li there is a term 
Ximfa^, which is of order m'jj/m'^ for i = 4,5,6 and additionally suppressed for i = 1,2, 3. 
The main originality of the is, however, the presence of a term with in and with 
^i,a j^. Computing the divergent coefficient one finds 



u 



327r2 



dx / dy (/C-lnAf) 



(7.6) 



where the infinite constant /C is defined in Eq. ( p.l4|) . Now, /C is independent of i and cancels 
out when one sums over i. This happens in the case of because Xi = 0, and in the case 
of M*'" because (LaL)j)^r = (rarD/^r = 0; the last relation holds because the matrices are 



diagonal, which is a specific property of our model. In u*, Eq. ( [7. 61 ), and also in one may, 
therefore, discard the divergence. We may also substitute lnA|^ by the logarithm of the 
dimensionless quantity /AY This logarithm is large for i = 4,5,6 while for i = 1, 2, 3 
it is of order mjy/ {mj^mjp). The order of magnitude of J^i^i"^^ Xljl-^a)/^*!-^!)*!-^*'" is, 
therefore, m'jj/mj^!^ 

In summary, J2i,=i ^Zi, suppressed by m^^/m^, even in the presence of charged scalars. 



B. Graphs in which the Z attaches to neutrinos 

There are also contributions to in which the Z boson attaches to the neutrino line, 
thereby changing the neutrino eigenstate from Xj to Xi- The corresponding vertex is given by 
Eq. ( p.l5|) , with an extra factor 2 because of the self-conjugated character of the neutrinos. 
One obtains a contribution to M^, and = X]t=i Let us define 

= xm] + yrn^j + {1 - x - y) (m^ - xm^ - ym^) - xyq^ , (7.7) 

and similarly D'^j, which is identical to with mw substituted by ma- We then find 



w 

ij 



UIUl 



There is an infrared divergence when rrii = 0. Therefore, it is better to take rrii = mi as the 
subtraction point for the logarithm. 

^We remind the reader that we do not take into account the enhancements introduced by large 
logarithms. 
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Ur 



-^9 



1 fl-x 



+ r 



- (Ra),. {UlUl)^^ (^I),.7l] (In/^r, - /C + 1)} 



with 



^1 = mimjYlL , 

^2 = [2pi ■ P2 (1 - a;) (1 - y) + mjx (a; - 1) + mly {y - 1)] 7''7i 
+m^mr (1 - X - 7^7r 

+2m^ (1 - x) \p^x - (1 - y)] 7l + 2m^ (1 - y) [p^y -p^{l- x)] , 



(7. 



(7.9) 



(7.10) 



+mimr (1 - y) Y 
+miX {m^Y - 2^2) 
+mjm^ {1 — x)Y 



{ulU,)^^ (4)^.^7^ - (LJ^, (UlUl)^^ {Ri),^lR 



+mjy 

+m^mj. {1 — X — y)Y 
+m^x (1 — x) 



X {m,Y - 2p^) 
+mry{l-y) 



X 



{nirY - 2p?) 



(Ra),. (UlUl)^^ {Ri),^lR - (La),, (uIUl),, (4)^.^7, 



+xy 



(Ra),. {UlUl) {Ri) (2m^K7L + 2my,iR - 2p, ■ p,Yil) 



(7.11) 



- (LJ^^ [uIUl),, {Li)^^ (2mX7ii + 2m.p^7L - 2p, ■ p,Y1r) 

As before, fC does not depend on rrii and on rrij and cancels out when one sums Eq. ( [7.8|) 
over i and over j. The cancellations occurs in the second line as a consequence of Eq. ( p.6| ) 
and in the fourth line as a consequence of (|2.8| ) and the form of Ra in Eq. ( p.l8|) , which are 
properties of the general seesaw framework; in the case of the term [La),i{U\UL)ij{L'^a)jT^ 
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cancellation upon summation over i and over j also hinges upon the fact that the matrices 
Ta are diagonal, which is a property of our specific model. 

We are therefore free to subtract from \tlD^ in Eq. ( |7.8| ) its value when i = j = 1, i.e. 
its value when both Xi and Xj ^ire the lightest neutrino. We obtain In (D^/D^), which is 
zero for z = j = 1, of order m%/ (m|jm|n) when i G {1,2,3} and j G {2,3} or vice-versa, 
and large otherwise. We similarly substitute InDfj by ln(Dj" /D"^). 

It is now possible to evaluate the order of magnitude of each contribution to M^^. After 
tedious yet straightforward consideration, we conclude that all terms are suppressed by, at 
least, m?jj/m\. In most cases this suppression applies term by term; in some exceptional 
cases one must sum the contributions over the light neutrinos and apply Eq. ( |6.25D . This 
happens, for instance, with 

/I /* 1 — X "1 

dx dy-^, (7.12) 

which is of order rrfp^ for i,j G {1,2,3} but acquires a suppression mj^/m'j^ when i and j 
are both summed over the light neutrinos. 

In conclusion, the contributions to the decay ( |3.1| ) from both photon or Z exchange are 
all suppressed by mj^/mj^, even in the presence of extra scalar doublets. This suppression 
also ensures that the simpler processes yU^7 and Z r+yU~ are invisible in all feasible 

experiments. We next look to the contributions to the decay ( p.l|) from neutral-scalar 
exchange. 



VIII. T~ fi-S^* 
A. Self-energy graphs 



Similarly to what happens with the couplings to the photon and Z boson, there are two 
self-energy graphs for the coupling to the neutral scalar S^: 

M,i = ^u, (r,7L + tI^r) H^Pi)] , (8.1) 

M,2 = ^u^ H^iP2)] ^\ ^ [VblL + Tl^i^ Ur. (8.2) 
V 2 mf^ — \ J TT 

The difference relative to the case of the gauge bosons lies in the fact that couples 
differently to the charged leptons r and /i — in the first case with the Yukawa coupling (r^)^^, 
in the second case with ^h)^^- As a consequence, it is not possible to use Eqs. ( |6.10| )-( p.l2| ) 
and we must compute S(p) explicitly. Define 

Df = xm\ + (1 - x) (m^ - xp^) , (8.3) 
D1 = xm^ + (1 - x) {ml - xp^) . (8.4) 

The relevant integrals are then 
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/f = /'rfxln^g, (8.5) 







1 Df 



gt= I rfx(l-x)ln-^, (8.6) 







together with fj^ and analogously defined. One has 

- 5^ />. [{L^)^, (Ri)^^ 7L + (i?a),, (4),,. 7k] I . (8.7) 

a ) 

As in the previous section, in the original definitions of /f and there should be a divergence 
/C added to the logarithms. However, JC is i-independent and yields a null contribution to 
S(p) upon summation over i. This happens because of Eqs. (p.6|) and ( p.7|) and because 
in our model the matrices Ta and are diagonal; in the second line of Eq. ( p.7| ) we also 
need the first Eq. ( ^M for (LamRD^r = (r.t/imt/^A,)^, = (E^MI^A^)^, = 0. In the 



definitions of f"" and of gf we could, therefore, subtract from \nD^ its value for i = 1, i.e. 
for the lightest neutrino; this subtraction corresponds to an ^-independent subtraction from 
/f and from g1 and, therefore, to a null contribution to ^(p). 

It is clear from Eqs. (|8.1| ) and ( |8.2| ) that the relevant p"^ in S(p) is either or m^. When 
i = 1, 2, 3 one has <S <C m^; when i = 4, 5, 6 one has <C -C m^. It is clear that, 
for i = 1,2, 3, the logarithm of Df /Df is — neglecting -C — approximately proportional 
to {mf — mf) /m\ ~ ^'d/ (^/j^f)- The same happens with the logarithm of D^ / DY ■ 
Therefore, the contributions of the light neutrinos to S(p) are suppressed by m|)/ {mj^m'jp). 
On the other hand, some contributions of the heavy neutrinos remain unsuppressed. For 
i = 4, 5, 6 one has 

2 

77? ■ 

/r^ln^, (8.8) 
9t^U^^—,-A ■ (8.9) 



2 V m 



Matrix elements of Ul, as in Xj and in {La)^i{L'^a)ir, suppress some of the contributions of 
the heavy neutrinos. Overall one obtains the unsuppressed part 

-m^ [(La),, K),,7L + {Ra),, (4),,.7i?] } . (8.10) 

Note that the sum over a includes a contribution from the charged Goldstone boson. 
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FIGURES 



Pi / k+pi \ pi p2 
— ' I — * — I — »— 



Pl P2 '' ^^+P2, P2 
1 — » 1 — ' — ' — »— 

Xi [I' 



-<- . 

k 



Pi P2 



5Cj 



FIG. 1. The three diagrams for r —>■ fi S?* which have parts unsuppressed by 



m 



R ■ 



We conclude that some contributions of the heavy neutrinos to M^i and to Mh2 remain 
unsuppressed when niR — oo. Let us compute those contributions in detail. Using 

maUL)u = {Mh)eeiUn),,, (8.11) 

which is valid in our specific seesaw model for any i and i, we rewrite Eq. ( p.lO| ) as 



1 ^ 

167r2 5Z 

i=4 



R 



fii \ R 



Ui) In 
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(ra),, (M^)^^ (AJ^, 7L - (A:)^^ (Mnl^ (Kl^ 7r 



fj.fl 

6 



MM 



IGtt 



i=4 >- 



(AD.. (Aa 



(r^),, (M^)^^ (A,)^^ 7L - (A*,),, (Mnl^ (F*)^^ 7i?^ 



'MM 



MM 



(8.12) 



where /i is an arbitrary mass, inserted for dimensional reasons; the expression ( |8.12| ) is 
independent of it. Inserting this result into Eqs. (|8.1|) and ( ^.2|) , we obtain — see Fig. |I| for 
the relevant self-energy graphs — 



with 



— 



Uf, {AlJl + Ar^r) Ur , (8. 13) 



Aj 



(Efe)^^ - (Ffe), 



T 



+ 



171^- (F 



&Jrr - "^M 



(r.)J W)^^5^(F,)^^(A, 



(8.14) 



A 



R 



(Efo) - (F^)^ 



^ E (AD,, (A.).. - (M^),, J] (F,)^^ (A, 



+ 



rrir (F 



6;rr - "^M 



(r^)J (^i^uE(^D,,(n 



(8.15) 
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B. Graphs in which 5^ couples to charged bosons 

The necessary couplings are found in Appendices |A4| and |A 5|. Let us define 



Df = xml + yml, + {l - x - y) {ml - xml - yrr?^ - xyq^ 
Then, the graph in which couples to (incoming) S*^ and 5^ yields 

1 



(8.16) 



—I 



EEC 



1 pl—x 

aa'b I dx dy 



-rrii 



Ur ■ 



(8.17) 



The cases in which either S*^, or 5^, or both, are charged Goldstone bosons, are implicitly 
covered through Eqs. ( |2.26| ), (|2.27|) , and ( |A34| ) , (|A36|) . The graph in which 5*^ couples to 
W~ and gives 



32^2772 



1 rl-x 



^ ""m \ nwa [(^ + 1) ^rlR + (y - 2) m^7i] 



M 1 jjWo 

[{2x + y) m^rrir'jR + x [1 + x) mlji 



D 



Wa 



+y {y - 2) ml^L + 2x{y- 2) pi ■ p^Il] - 2 (LJ^, In 
Similarly, the graph in which couples to and contributes 



haW 



327277^ 



1 rl-x 



Jo 



X U, 



rrii 



D 



aW 



[{y + 1) m^7L + (x - 2) nirjR] 



(Li) 



D 



^ [(2?/ + x) rrif.m-r'^L + x{x-2) ml'^R 



Finally, there is a graph with the 5*° attaching to two bosons, yielding 



Ml 



bWW 



ig^2 ^1^^ 2^^^ — — 



(8.18) 



+y{l + y) ml^R + 2y{x-2)p,- p^Ir] - 2 {L\) 7^ In . (8.19) 



(8.20) 
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Any part of the integrals which does not depend on rrii ends up, upon summation over i, 
giving a zero contribution to M^. One may therefore subtract from each integral its value 
when rrii = nii. In Eqs. ( |8.18| ) and ( p.l9| ) we have already performed that subtraction 
in the logarithmic terms. The infinities which occurred together with the logarithms in 
those expressions have been dropped using [LaU'D^r = iXa)iJLT = in Eq. ( |8.18| ) and the 
Hermitian-conjugate relation in Eq. ( ^.19|) . The amplitudes (|8.17|) and (|8.2CI| ) have only 
finite integrals. 

Consider for instance Mi^aa' in Eq- ( |8.17|) . The term with Ra'Ri is proportional to 
{Ur) ^i{U\^iT-. For i = 1, 2, 3 it is suppressed by m\,/m\ from the matrix Ur, it is addition- 
ally suppressed by Caa'b and the integral, giving together a factor mj:,/mp. For z = 4, 5, 6 
the suppression factor is mpmri/rn?^. 

The term with LaiL\ has {UL)i_d{Ul^)iT- For i = 1, 2, 3 we neglect the neutrino masses rrii 
in the integral and use Eq. ( |6.25|) to obtain the same order of magnitude as for the term 
Ra'Ri- For z = 4, 5, 6 we have suppressions mjj/m^ from Ul and mprriD/mj^ from the rest 
of the term. 

The terms with Ra'L\ and La'R^ have suppressions niD/'mR from the mixing. For i = 
1,2,3 there is also rrii and the integral, yielding together suppressions m'jj/ {mRrriF)- For 
i = 4, 5, 6 mj ~ rriR and the integral is of order rn^ ^ and we correspondingly obtain mp/rriR 
in addition to the mixing suppression. 

In conclusion, all contributions to M^aa' are suppressed by, at least, rripmu/m'j^. In 
the same way, one finds that M^wa and M^aw are suppressed by m\,/m\ while M^ww is 
suppressed by m|,/ {m\rriF)- 



C. Graphs in which 5"^ couples to neutrinos 

We remind the reader of the quantity defined in Eq. ( [7.7] ), and of the analogous 
quantity D^^ We find, for the contribution to in which 5'° couples to neutrinos, the 
result 



M, 



Here, 



167r 



^0 



(8.21) 
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7W 



/ iUL)^i (Ul^^^ (1 - x) m^m, (fI^^^ 7l + (1 " v) ^i^r {F,)^^ 7^ 
-ymjrrir (^F^^ .iR- xm^rrn {Fb)-j 7l , 



(8.22) 



Z"- = —rriim 



+m,imr (1 - y) 



—rriim^x 



(Ra),. [Fn {Ri)^^^p + (La)^, {F,\^ {Li)^^^R 
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+m^mj (1 — x) 



-mrmjU 



+ [x (1 — x) + ?/ (1 — y) ml — xy (m^ + ml — g^)] 



and 



(8.23) 



(8.24) 



In Eq. ( p. 211 ) we have already dropped the infinity occurring together with the logarithm, 
because, using all the previous arguments for the cancellation of terms independent of m^, 
we find Y.i,j = 0- 

It is tedious but straightforward to check that all the terms in Z'^ and in Zf^ end up 
suppressed by 1/"^!; or by higher powers of m,~j^. The same does not happen, however, 
with the terms in (which include a contribution from the charged Goldstone boson for 
a = aw). Let us write the latter in more detail, using Eqs. ( p. 181 ), ( p.l9| ), ( |2.22| ), and the 
fact that the matrices Fq, A^, and Af, are diagonal: 



2^:" = (ra),, (Aa)^^ {ulAlUn + U^AlUl)^^ (uj,) 

+ (K),, iK)rr (Un),. {U^A^U, + f/jA,f/*)^ (f/t) 



(8.25) 



Due to Eqs. ( [4.1[ ), terms unsuppressed by mn/mR may arise only from the first term in 
the right-hand-side of Eq. ( ^.25|) , when z = 1, 2, 3 and j = 4, 5, 6, and from the third term 
in the right-hand-side of Eq. ( ^.25] ), when z = 4, 5, 6 and j = 1,2,3. In the first case the 
integral of InD^- is practically ^-independent, in the second case it is almost j-independent. 
One obtains the following contribution to M54, with the corresponding Feynman diagram 
depicted in Fig. |l|: 



64 



16v^7r^ 



E 

j=4 



R)l,i \ ^R 



UU In^ 



+ (A,),,5^(A*)^^(F* 



kJrrlR 



U, 



■7l 



(8.26) 



This expression is independent of the arbitrary mass /i. Mfe4 is not fully suppressed by 
powers of since M^^ becomes constant in the limit m/j — > 00. 
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D. Unsuppressed terms 

We thus conclude that the vertex — fi^ S^* does not vanish in the hmit mji oo. 
One has 



16V2' 



i=4 



2 ■ 

mf 



m?, — m? 



7l 



+ 



/I 



m?, — m'i 



IR > Ut , 



(8.27) 



with Al and given in Eqs. ( p.l4|) and ( |8.15| ). Note that is independent of the 
arbitrary mass parameter /i. 

It is interesting to observe that Mb is suppressed when there is only one scalar doublet. 
Indeed, in that case there is only one physical scalar, the Higgs boson, which has 6=1. 
Moreover, there is only one matrix T = *) Ml and, smce [\/2/v) Md = A, we find 



A, = ^ {ml - ml) (MB),^. (F)^,^ (A), 
An = ^ {ml-ml)iMnUiA*)^^iT* 

Thus Eq. ( p.27| ) vanishes in this simple case. 



Mfe°° = . 



(8.28) 



IX. BOX DIAGRAMS 



There are four classes of box diagrams, as depicted in Fig. ||. In that figure, either S^, 
or S^,, or both, may be substituted by W^. 

Let us first consider the box diagrams in which the fermion line starting in the incoming 
T~ ends in the outgoing Those are the diagrams denoted in Fig. g "box 1". The 
box diagram with two charged scalars and S^, gives the following contribution to the 
decay (§3): 



Maa' 



167r2 



1 rl — x 

dx / dy 
Jo 



l-x-y 



i,j a,a' 

The box diagrams with one 5*^ and one give 

_,2 ^ s. /"l pl—x fl—x—y 



A aa 



Tjaa 



M, 



aW 



Wa 



327r' 



327^' 



dy 







dzj2 



AaW 



r)a\ 



aW 



9 piaW 



A 



Wa 
ij 



B. 



Wa 



(ATI' 



9 p\Wa 



(9.1) 



(9.2) 
(9.3) 
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<- 
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s+ 
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P4 



5Ci 



P2 



1^ 
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Pi 



P3 



Xi 
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Pi-k 
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k-p3 

box 4 



5C, 



P2 
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I I ( i . 2. The four types of box diagrams for r n e^e which occur in our model. 
Finally, the box diagram with two yields 



- ^ ^^^^Mi V"^A. V"^ 

/•I pl—x pl—x—y 

X dx dy dz 
Jo Jo Jo 



A 



WW 



B. 



WW 



^D^wy 2D 



WW 



(9.4) 



In Eqs. O-dl)' 

D^f = xml + yml, + zm^ + {1 — x — y — z)[ 



2 2 2\ 

- xp^ - yp2 - zr ) 



-xyq^ - xzpl - yzpl ; 



(9.5) 



D"^, D^'^', and D^^ are obtained from D^^ by substituting either rua', or m^, or both, 
by rriw In order to write down the expressions for A""', . . . , one must define 



One then has 



Ada _ - 
^ij — 



P'" = p° + ;Z - 1) + + p2 (1 - Z) 



(9.6) 
(9.7) 
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A 



Wa 
ij 



-u,. 



rrin 



AW 



WW 



(9.8) 

(9.9) 

(9.10) 
(9.11) 



One obtains 3^°"' by discarding from all the terms which contain neither /P nor Z^', 
and then substituting, in the remaining terms, P by 7^ and /" by 70. Applying the same 
algorithm, one obtains Bf^ from B^"" from AY,"", and Sjf ^ from A,^^. 

Next we consider the box diagrams in which the fermion line starting in the incoming r~ 
ends in the outgoing e~ . Those are the diagrams called "box 2" in Fig. ^ The corresponding 
amplitudes may be obtained from Eqs. (|9.1| )- (|9.11| ) by using, instead of Eqs. ( pl5|) -( 



2, 2, 2, /I \ ( 2 2 2 2\ 

xm^ + ym^, + zrrij + [l — x — y — z) [m^ — xpi — yp^ — zq ) 



-xyr 



xzpl - yzpl 



u. 



p-=p'^ ^x + z)-p'^z + pty, 
P'-=p- (x + z-l)+pUl-z)+pty, 

and by making, in Eqs. ( |9.8|) -( pnD , the substitutions 
{La,a> ) Ra,a') {La,a' , Ra,a') Finally, ouc must chaugc the overall sign of the amplitudes, 
i.e. insert a minus sign in front of Eqs. ( p.l|) -( p74D , due to the interchange of two fermions in 
the final state. In Eqs. (^]2|)-(^]4|) one must also interchange {Ul)^^ — > (^i)ej i^L)ej ~^ 
{Ul] 



(9.12) 
(9.13) 
(9.14) 

Me, {La,Ra)ei ~^ {La,Ra)nj^ and 



Another type of diagrams are the box diagrams denoted "box 3" in Fig. ^ Those 
diagrams arise due to the Majorana character of the neutrinos, and they must be computed 
using specific Feynman rules for Majorana fields — see, for instance, Ref. fl^. One obtains, 
analogously to Eqs. ( |9.1| )-(P^), the following contributions to the decay (|3.1|): 
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l-x-y 



(9.15) 




(9.16) 



(9.17) 
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X dx dy dz 
Jo Jo Jo 



^ww 



B 



WW 



f)Ww\ 2D 

ij J 



WW 



(9.18) 



Here, 



7-jaa 



xml + yml, + zw?, + {1 — x — y — z) 



-x{l- X - z)pI- y{x + z)pI- y{l- y - z)pI- x{y + z)pI 

+xyq^ + [xy — z {1 — x — y — z)]r'^ . (9.19) 

Notice that Mww is proportional to rrii and to mj, indicating that it vanishes in both hmits 
rrii —>■ and mj 0. This is an instance of Kayser's "practical Majorana-Dirac confusion 
theorem" [|l^. Defining 



one may write 



A aa 



P'- = p^ (a; + z-l)+p'^y + p2{l-z) 



(4). (LD JlR+Ul), {R\)JlL-mARi), {Ri),^7L 



(9.20) 
(9.21) 
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~\aW 



(i?J,,/" + (Lj^,m, 



Cut 



(9.22) 



(9.23) 



A, 



Wa „.T^-1 



viC- 



A 



WW 



(9.24) 
(9.25) 



One obtains Bff from iff'— and similarly Bf^ from iff, M"" from A^", and from 
AjJ^ — by deleting all the terms which contain neither /' nor and then substituting, in 
the remaining terms, f by 7^ and /" by 'jg. 

Finally, there are the box diagrams of the type denoted in Fig. ^ "box 4". The corre- 
sponding amplitudes may be obtained from Eqs. ( |9.15 )-( |9.25| ) in the following way. Firstly, 
one must change the overall sign of Eqs. ( |9.15|) -( ^.18|) and perform in them the inter- 
change {Ul)^- {Ul)^-. Secondly, in Eqs. ( p.l9[ )-( p.21[ ) one must interchange P2 ^ P4, 
vi ^ P45 cind 



^ r 



'. Lastly, in Eqs. (|9.22|) -( p.25|) one must interchange 



u, 



^ _ ^ Me and 

{La^a'iRa,a')f.j ^ {La^ , Ra,a') ■ 

One may analyze the mji dependence of the box amplitudes given above by using the 
skills developed in the previous sections. It is easily concluded that all those amplitudes are 
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suppressed by at least one factor m^^. In order to reach this conclusion, it is necessary to 
use Eq. ( |6.25|) , when the box diagrams are either of type "box 1" or "box 2" and both Xi 
and Xj ^^^6 light neutrinos. 



X. THE ASYMPTOTIC LIMIT rriR ^ oo 

In our model there is a scale m^j which is much higher than the other two scales, m£, and 
TUp (see Section fV\j . In this section we want to study the asymptotic limit m^j — oo of our 
model. The simplest way — the one which we have in mind in the following — of increasing 
the scale mji is by multiplying the mass matrix Mr in Eq. ( |2.3| ) by a (dimensionless) factor 
which becomes much larger than one. 

When the decoupling theorem |[1^J2^ applies straightforwardly, one may simply delete 



the heavy fields from the Lagrangian in order to obtain the low-energy theory. In our case 
this does not work because, if we remove the fields from the Lagrangian, we also delete 
any trace of the flavor-changing neutral interactions. This is at odds with the explicit one- 
loop calculation of the vertex i i'S^* , for i ^ i', since that vertex does not vanish in the 
limit mn — > oo (see Fig. [l| for the Feynman diagrams with non- vanishing contributions). 
On the other hand, according to Ref. the limit oo must yield a sensible theory. 

Evidently, the only theory which can emerge from our model in that limit is the multi- 
Higgs-doublet SM, containing flavor- changing neutral Higgs interactions suppressed by loop 
factors IGvr^ and by small couplings. Thus, we expect a Yukawa Lagrangian of the form 



= - E ( ^7 ' ) H.c. (10.1) 

In order to demonstrate that this Lagrangian indeed emerges, we consider the following 
couplings: 

1. iL^i'^S'*, 

2. UL ^ iRS+*, 

3. ul ^ 

For the sake of brevity, let us denote the corresponding tree- level vertices by Vq, V+, and 
Vw, respectively. The one-loop contributions to these vertices fall into three categories: 

I. The contribution has a non- vanishing limit mR oo. 

11. In the limit oo the contribution vanishes. 

III. The contribution is independent of and is present also when the fields z/^ are 
removed from the Lagrangian. 

Our strategy is the following. We identify all contributions of Category I, since only they 
are relevant for obtaining the limit mn — >■ oo of our model. We then show that there are 
three contributions of Category I to Vw, but they cancel out except for a part which may be 
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viewed as a unitary transformation on the vector of flavor fields vl- The transformed field 
vector is denoted by ul- After these steps, we see that the contributions of Category I to 
the vertices Vq and V+ are identical, provided that the neutrino field is used in V+. This 
demonstrates that in the limit — oo the Yukawa Lagrangian of Eq. (|10.1|) emerges. After 
the one-loop corrections, as discussed above, we remove the fields from the Lagrangian. 
Then all the traces left by those heavy fields are contained in the off-diagonal elements of 
Fj, which have arisen from the one-loop corrections of Category I. 

We now pursue the strategy outlined in the previous paragraph. First consider the vertex 
and self-energy corrections for the above couplings: 

H : Tfc + AFfc (H, S^) + ^ [4 {S^)]^ F, + ^ Ffczi {S^) , (10.2) 

V+:V, + AF, S') + \ [4 [S^)]^ F, + ^ V,zl [S') , (10.3) 

VW : ^ { 1 + F (W, 5) + i [4 {S^)]^ + ^ 4 {S') } . (10.4) 

For the definition of the wave-function renormalization matrices 4' 4' ^^"^ ^P" 
pendix It is proven by a tedious checking of all one-loop graphs that the corrections to 
Vo, V+, and Vw not included in Eqs. ( |10.2|) , (|10.3|) , and ( |10.4|) , respectively, are either of 



Category II or Category III. In the expressions ( |10.3| ) and ( |10.4| ), only one- loop corrections 
in which the external neutrino legs correspond to light neutrinos are interesting. Those 
one- loop corrections are obtained by performing the calculations with the fields Xii but the 
result has to be multiplied appropriately by Ul] in order to obtain AFfc(V^+, S^) a multipli- 
cation from the right by U\ is necessary; the 6x6 matrix part of the neutrino self-energy 
associated with the Dirac structure has to be multiplied by Ul from the left and 
from the right in order to arrive at the 3x3 matrix we need; and so on. With this pro- 
cedure we have in mind that, in the limit m^j 00, the light neutrinos become massless, 
and we are allowed to work with the fields vli which are members of the left-handed lepton 
doublets in the Lagrangian before spontaneous symmetry breaking. In the following, our 
results for the quantities appearing in Eqs. ( p.0.2| ), (|10.3|) , and (|10.4| ) will only be given for 



rriLi 00, since this is the limit we aim at. In Eq. (|10.2|) the vertex correction AFfc(Vo, S^) 
is obtained by exchanging all charged scalars, including the charged Goldstone boson G^. 
In Eq. (|10.3| ) the vertex correction AFfc(V+,S'°) stems from exchanging all neutral scalars, 
including the neutral Goldstone boson G^. The vertex correction F(VV, S) originates in the 
coupling of the W boson to a charged and a neutral scalar; thus, the loop has S^, S^, and 
Xi- Concerning the wave-function renormalization matrices 4(*S'^)) 4('^^)' -^l(*S'°), in 
parentheses we have indicated the scalar exchange they come from in the self-energies. 
In Eqs. ( |10.2| ), ( |10.3| ), and ( |10.4D , all the one-loop contributions which occur contain two 



A and one F Yukawa-coupling matrices. The vertices Vq and V+ also receive contributions 
from the wave-function renormalization matrix of the Higgs doublets; these contributions fall 
into all three categories. One can show that those of Category I induce the same corrections 
to Ffc at the vertices Vq and V+. 

Before we list the one- loop results for the quantities appearing in Eqs. ( |10.2| ), ( |10.3|) , and 
( |10.4| ), we want to introduce some useful notation. First we note that 
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lim Ur = (0, W) , lim mi,2,3 = 



(10.5) 



where is a 3 x 3 unitary matrix. We define a useful 3x3 Hermitian matrix by 



2 4 2 



10.6) 



wliere tlie divergent constant /C has been defined in Eq. ( ^.14[ ) and 

rh = diag {1714, m^, ttiq) . 
Then, it is easy to write down the self-energy 

for the neutrinos, and 

{s^;p) = Mil + b^l + b^^r 

for the charged leptons, with B given by 



B 



167r2 



\ M^^{K, + l)-Wm\nm'^W^ 



A, 



:io.7) 



:io.8) 



(10.9) 



(10.10) 



The superscript 00 reminds us that we have taken the limit ttir — 00. In the same limit, 
the vertex corrections are given by 

ATk {Vo, S^) = ATk {V+, 5°) = ^ 5^ r,At (/C + 1 - Wlnm^W^) A, (10.11) 



and 



T{Vw,S) = -A. 



(10.12) 



The off-diagonal elements of Arfc(Vo, S^) are given by the vertex correction computed in 
Section |V111| , in the limit rriR 00; the same holds for Tj'^{S^;p). In order to calculate 
the self-energy ( |10.8|) one has to take into account the Majorana nature of the fields xf- the 
neutrino self-energy derives from a propagator at second order in perturbation theory, and 
there are 2^ = 8 possibilities to attach the external legs to neutrino fields in the two 
Yukawa Lagrangians; therefore there is a combinatorial factor 2^/2! = 4. Furthermore, the 
terms Bl,r drop out of the neutrino self-energy because of [|1^] 

J2bkbk' = 0. (10.13) 



We have also used the relations (see Appendix A.l) 
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Kbk' = 2 ^ alak' = 25kk> 



(10.14) 



which have enabled us to sum, in Eqs. ( |10.6|) , (|10.10 ) and ( 10.11 ) over the index j. 

Since in the hmit m^j — > oo the neutrinos are massless, the procedure for the renormahza- 
tion of the self-energy laid out in Appendix B.l is not applicable. However, it is reasonable, 
in view of Eq. ( |10.8[ ) , to define 



(10.15) 



The determination of and of follows the on-shell prescription in Appendix B.l. 

In the expression for the one-loop correction of IV, Eq- ( |10.4D , the only matrix which 
is possibly non-Hermitian is z\^[S^). In any case, we may decompose that matrix into a 
Hermitian and an ant i- Hermit ian part: 



with [zi)^ = ± 



;io.i6) 



We notice that in the T,^{S^]p) of Eq. (|10.9| ) there is no term A^^^. Using the equations in 
Appendix B.l for calculating the fermion wave-function renormalization matrices, this fact 
simplifies considerably the expression for z\^ ^. Since A is Hermitian, we obtain 

Re[4(5±)-4(^±)]_ = A„„, Im[4(^±)-4(^±)]_ = -— ImS„„, (10.17) 



and 



ReK(^^)+4(^^)].„ = ^o 
Using the possibility to choose the convention 



(10.18) 



(10.19) 



as explained at the end of Appendix [B 1| , we obtain [^^(S'^)]^^ = (-2+)aa = ^aa and {z[_)aa = 
0. For a 7^ /5 we apply again the procedure of Appendix B.l and obtain {z^)f3a = A^a- In 
summary, with the convention of Eq. ( |10.19|) we find 



A 



and, with Eqs. ([TaT2|) and ( |Tal5| ), 



1 



r(vw,^) + -(4 + ^L) =0. 



Thus, Eq. ( p.0.4|) is given by 



(10.20) 



(10.21) 



;i0.22) 



Since the light-neutrino masses vanish when m/j ^ oo, we are allowed to rotate unitarily. 
An infinitesimal unitary rotation is given by 
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ul = (1 + n) ul , with = -n. (10.23) 



With the choice |2T| 

1 



1^=^/, (10.24) 



the one-loop-corrected W vertex of Eq. ( |10.4| ) reduces to the trivial form g / ^/2, if conceived 
as pertaining to ul- 

In terms of the new field ul, the vertex V+ is given by 

r, + Ar, (y+, 5°) + i [4 {s^)]^ + i r, [zi (5°) + z'_] . (10.25) 

Since 

{S^) +zi = A + zt = zi {S^) , (10.26) 



and because of Eq. ( p,0.11| ), the expression ( |10.25| ) for the charged- scalar vertex is identical 



with the expression (|10.2|) for the neutral-scalar vertex. Therefore, in the ul basis for the 
neutrino fields, the one-loop corrections — associated with — to the couplings of types 
^r'^lS^ and IrVlS' are identical. This is crucial for writing the theory after decoupling of 
the heavy fields as a multi-Higgs-doublet Standard Model with the Yukawa Lagrangian of 
Eq. (|10.1|) . The coupling matrices are then given by Eq. (|10.2| ). The infinities introduced 
by the one-loop contributions of Category I are all in the diagonal — see Eq. ( [10.11| ) — and 
they may, therefore, be absorbed by redefining the diagonal matrices as renormalized 
coupling matrices. This concludes our argument, valid at least at the one- loop level, that 
in the asymptotic limit niR oo one obtains the Lagrangian of Eq. ( |10.1| ) and our model 
approaches a multi-scalar-doublet SM with suppressed off-diagonal couplings in Vj. 

Several remarks are in order. First, we stress that the contributions to 4 from ex- 
change, and those to 2;^ from exchange, are flavor-diagonal and Hermitian and, therefore, 
they do not have anti-Hermitian components which would interfere with the arguments pre- 
sented above. This must be so because these contributions belong to Category III. Second, 
though in the limit of infinite right-handed scale we were able to show — taking into account 
a rotation of vl — that the one-loop contributions of the right-handed neutrino singlets are 
the same for couplings of the types IrIlS^ and IrUlS' , this does not happen with contri- 
butions of Category III, due to the different mass effects of charged leptons and of massless 
neutrinos; thus, the fully one-loop-corrected coupling matrices do receive different finite 
parts in the iRiiS^ and d-RVhS^ couplings because of contributions of Category III, an effect 
which is to be expected in the multi-scalar-doublet SM. Of course, the infinite corrections 
to couplings of both types are the same, which allows for a consistent renormalization pro- 
cedure. Third, because of the cancellation in Eq. (|10.21|) , the infinities at the vertex Vy/ 



stemming from scalar corrections of Category I also cancel. This is necessary for consistency, 
since it would be impossible to absorb those infinities into the gauge coupling constant. 
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XI. DECAY RATES 



Experimental bounds on lepton-flavor-changing processes are found in Ref . |^ . For the 
r the bounds on the branching ratios are of order 10~^, but for n decays they are five to six 
orders of magnitude better. Of course, all our formulas are easily adapted to muon decays. 

In this section we shall always neglect the masses of all final-state fermions. Pursuing 
the philosophy laid out in Section |^ we shall also assume that Yukawa couplings are of 
order Y ~ rriD/mp, cf. Eq. ( ^.2| ). 

First we consider the process 7- Its matrix element may be written 

A {t~ fi-^) =ee*pU^iaP^qx {AI^l + ^^7h) , (11-1) 
where e is the polarization vector of the photon. Then the decay rate is 

r(r-^^-7) = r,^g|(Kr + l^jr), (11.2) 

where a is the electromagnetic fine-structure constant, Gi;' is the Fermi constant, and 

r.^r(r-^e-p,..) = ^- (11-3) 



The amplitudes in Eq. ( |1 1 . 1| ) are given in our model by Eqs. ( |6.17D and ( |6.18| ): 



KR = Y.^L,m{q' = 0). (11.4) 

i 

In Subsection [VI B| we have made for those amplitudes the order-of-magnitude estimate 
AJ^ ~ Y^m£)/(167r^m|j). Thus, the branching ratio should be of order aV^/ {TiG'^^mj^), 
where we have used rrir ~ niD. With Gp ~ 10~^ GeV~^ and m/j ~ 10^° GeV, we see that 
the branching ratio would be ~ 10"'^^ even if we allowed Y to be of order 1. 

Next we consider the lepton-flavor-changing Z decay Z T^fi~. Its decay amplitude 
may be written 

A{Z^ r+/i-) = ^ ePu^^, {Af^L + ^Itr) Vr , (11.5) 



where e now denotes the Z polarization vector. The decay rate is then 

r(2--v)=i^(i^!r+i^in- (11-6) 



The amplitudes Afj^ can be read off from J2t=i ^Zl^ computed in Section |VI1] . As they are 
of order mj^/mj^, once again one finds a ridiculously small decay rate. 

We now discuss neutral-scalar decay, which is not suppressed by inverse powers of ttir 
whenever uh > 2. The matrix element is 

A ^ r+/i-) = u, {Ai^L + A%jn) Vr • (11.7) 
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Then the decay rate is 



r - 

In our model we may identify 



levr 



[11.8) 



11.9) 



with Mj^ 



given by Eq. ( ^.271 ) without the spinors. This is not suppressed by any power 
but it contains three Yukawa couphngs. Thus, the decay rate should in general be 
very small due to a factor y®/(167r^)^. 



of m. 



In any case this is not very interesting, since no 
fundamental neutral scalar as been observed up to now. 

Finally we consider our model process r~ fi^e^e~ . Rather general formulas for the 
decay rate can be found in Refs. ||TT| , p!2| . Comparing our result of Section [V^ for the photon 
sub-process with the results in Ref . , we conclude that this contribution has a branching 
ratio of order [a^/(G|.m|j)](mD/'"^F)^/167r^. In the Z sub-process (Section |V11| ), exchange 
dominates over the exchange of charged scalars. Thus we estimate for the branching ratio 
of this sub-process [afij/{Gpmj^)]{mD/rnz)* (of course, mz ~ mp), where aw = g^/i^An) is 
the weak fine-structure constant. A similar suppression is found for the box sub-process of 
Section Therefore we may safely neglect all those contributions, including interference 
terms, and concentrate only on the neutral-scalar sub-process computed in Section |VI11 . 
The amplitude is then written 



A{t ^ /i e+e ) = ^ (Afe7L + p^.'^r) Ur (1^7^ + Tl'-fR) 



(11.10) 



which is modeled according to Eq. (|5.3| ). Thus, we are making the identifications = (Th)ee 
and, assuming more than one scalar doublet. 



Uf, {KlL + PhlR) Ur 



V2 



mt 



;ii.ii) 



where we have used the approximation ^ m^. From Eq. ( |11.10| ) one obtains 

r (r^ fi^e^e^) 1 
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;ii.i2) 



Both pb and A;, are of the form F^/(167r^m^), where mn is a typical neutral-scalar mass, 
which should be of order rriF- Therefore, the order of magnitude of the part of the branching 
ratio which is not suppressed by powers of rriR is F^/(167r^Girm|.)^. With IGn'^ Gpm'jp ~ 
10 and Y ~ 10~^ as reasonable numerical values, this branching ratio is quite small, of 
order 10"^*^, yet it is much larger than the photon, Z, and box contributions, which are all 
suppressed by m]^^. Due to the dependence on Y^, we can easily achieve a branching ratio 
much larger than 10~^^ by moderately increasing Y. However, in this case it is certainly not 
possible to allow for F ~ 1. In particular, if we apply the present estimate to the branching 
ratio of p~ — > e~e~^e~ , we rather find Y^ < 10~^°. 
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XII. CONCLUSIONS 



In this paper we have computed the amphtude of the lepton-flavor-violating decay r~ — *■ 
fi~e~^e~ in the context of the seesaw model with an arbitrary number uh of Higgs doublets, 
but with the assumption that the (tree-level) Yukawa couplings conserve lepton flavor. Our 
calculation is easily adapted to other lepton-fiavor-violating decays of the same type, e.g. 
H~ — >■ e~e~^e~; moreover, the parts of the calculation in which r~ fi~e~^e~ proceeds via an 
intermediate photon or an intermediate Z boson are also applicable to lepton-flavor-changing 
decays like /i^7 and Z — >• r^n^ , respectively. As a function of the right-handed scale 

TTiR we have found the following behavior of the decay amplitudes: 

A{t- fi--f) OC 1/771% , (12.1) 

AiZ T+fi-) OC 1/m^ , (12.2) 

, , 11/ m% for 7iH = 1 , / s 

A{r- ^ /.-e+e-) oc <^ ^ f ^ ' (12.3) 

I constant for 71h > 1 • 

The partial amplitudes for ji^e^e^ which behave like rn^ are, in general, suppressed 

by at least a factor mj:)7np/m\ ~ 10^^^; here we have assumed that the elements of the Dirac 
neutrino mass matrix Ma are of order m/) ~ 1 GeV, the Fermi scale is of order 7nF ~ 100 
GeV, and m^j ~ lO^'' GeV.0 The exception to this suppression occurs precisely when there 
is more than one scalar doublet, for then there are contributions to r~ — /i~S'°*, where 5*° 
is a neutral scalar, which are not suppressed by inverse powers of m^j. These unsuppressed 
contributions originate from exchange of the charged scalars S*^ ; the corresponding Feynman 
diagrams are shown in Fig. ^ We have found that the vertices yU~7* and r~ fi~Z* 

(and therefore also the decay amplitudes for r~ fi~'~f and Z —>■ r^fi'), which had already 
been computed before by other authors in the case of only one scalar doublet, are, even in the 
case of many doublets, suppressed by m|)/m|. ~ 10~^°. Furthermore, also the box diagrams 
behave like 1/771"^^ for an arbitrary tih-^ Therefore, with the above numerical assumptions, 
all classes of diagrams, except those with intermediate neutral scalars in the case of more 
than one Higgs doublet, are basically irrelevant for the process (|3.1| ). 

The unsuppressed contribution to the vertex r~ — fi~S^* is found in Eq. (|8.27|) , with 
the quantities and Aji given by Eqs. ( p.l4| ) and ( |8.15| ), respectively. That unsuppressed 



contribution has three Yukawa couplings and a denominator IGvr^; if we assume each Yukawa 
coupling to be of order 7710/771F ~ 10^^, we obtain an overall order of magnitude 10~*. This 



We want to stress that our philosophy is different from the one employed, e.g. in Refs. |TI,13|. 
We assume that m^j is of the order of or TUr, thus ttir is very large in order to implement the 



seesaw mechanism. In Refs. pl|Jl3| it is assumed that ttir ~ 1 TeV, or smaller, in order to obtain 
effects in lepton-flavor-violating processes. 

^Note that the processes t~ — > ^^/i^e+ and e^e^/^+ can, at one-loop level, proceed only 

via box diagrams [p^. 
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is certainly very small, but anyway much larger than anything which is suppressed by m'^. 
Moreover, we may allow for Yukawa couplings larger than 10~^. 

Estimating the branching ratio for the non-suppressed decay fi^e^e^, we obtain 

/ (IGTT^Gi^ml^)^ ~ lO""*^^, where we have taken uih ~ for a typical neutral-scalar 
mass, Y ~ mn/mp for a typical Yukawa coupling, and with the same assumptions as 
before for the scales irtp, and tjir. This branching ratio, valid for uh > 1, is not 
suppressed by the large scale rriR, but it is suppressed instead by the eighth power of a 
typical Yukawa coupling Y. We want to stress, however, that the estimate 10~^^ for the 
branching ratio of fi~e~^e~ is rather crude; having in mind the remarks at the end 

of the previous paragraph, it could easily be several orders of magnitude larger due to its 
dependence on Y^. On the other hand, the branching ratio for ^ yU~7 is suppressed by 
aYy{TrGlm%) ~ 10-^°. 

We have studied in detail the non-decoupling in the Higgs sector for m/j — > 00. In 
this limit our model approaches a multi-Higgs-doublet Standard Model with lepton-flavor 
non-diagonal Yukawa couplings, the off-diagonal couplings being, however, suppressed. The 
situation here reminds somehow the Minimal Supersymmetric Standard Model, which has 
two Higgs doublets, where non-decoupling in the Higgs sector has been found when the 
SUSY scale is made much larger than the Fermi scale but the masses of all scalars are kept 



of order mp 23 



The model discussed here with lepton-flavor-diagonal Yukawa couplings was put forward 
in Ref. as a framework for imposing large or maximal neutrino mixing. Here we have 
shown in a detailed way that, despite the soft breaking of the lepton numbers at the very 
large scale mji, the branching ratios of lepton-flavor-violating processes remain very small. 
We have identified the class of processes whose vertices are not suppressed by m]^. In 
this class the most promising example for future experiments is /i~ — > e~e~^e~] this non- 
suppression requires more than one Higgs doublet. Those vertices suppressed by lead 
to branching ratios far beyond present or future experimental limits. Thus we have a viable 
model with the interesting feature that neutrino mixing has its origin at the ultra-high scale 
rriR — the order of magnitude of the masses of the right-handed neutrino singlets, which is at 
the same time, via the seesaw mechanism, responsible for the smallness of the light- neutrino 
masses. 
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APPENDIX A: FORMALISM FOR THE SCALAR SECTOR 



1. 5± and 

We allow for an arbitrary number nn of scalar doublets 0^ {k = 1,2,..., tih) and use 
the notation 

yk J V 2 

We then write 

^1 = ^ + ^1', hence (0|^°'|0> = 0. (A2) 
The quadratic terms in the scalar potential are written 

(A3) 

The matrix A^^ is complex and Hermitian, while the matrices A and B are real and sym- 
metric; C is real but otherwise arbitrary. They are all uh x uh matrices. The eigenvalue 
equations are 

M'la = mla, (A4) 

°\lmbj \C^Bj\lmbJ ''\lmbj ^ ^ 

where a and b = Reb + ilmb are tih x 1 complex vectors. The orthonormality equations for 
the eigenvectors are 

^ala[ = a'^a' = 6aa' , ^{RebkReb[ + lmbklmb[) =Re{b'^b') = 6w , (A6) 

k k 

a^afc/ = Re 6fc Re 6fc' = Im &fc Im 6fc' = 5fcA;' , and Re 6^ Imfe^/ = . (A7) 

abb b 

The physical (mass-eigenstate) charged scalars and the physical neutral scalars 5*^ are 
given by |T|1 

S: = J2^lvt and S',=V2j2^e{bld'), (A8) 

k k 

respectively; or, equivalently. 



(A9) 



{vk + j:bhs'b)/V2 

Indeed, one then has 



2 

b 
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2. The mass matrices of the scalars 



The scalar potential is 

v = Y. ^''Ah + E ^^^^^ ('^^'^0 ('^I'^O ' (All) 

i,j i,j,k,l 

with 

= (/^ji) ; -^ijfci = Afciij , and Xijki = X*uk ■ (A12) 
The mass matrix of the charged scalars is obtained from the potential in Eq. ( [All[ ): 

Ml = fi^ + A, with Aij = J2 ^^ikl vlvi . (A13) 

k,l 

The matrix A is Hermitian. We proceed to find out the terms in V which are quadratic in 
the neutral scalars. We need to define two more uh x hh matrices, K which is symmetric 
(but, in general, not Hermitian) and K' which is Hermitian: 

Kik = ^ Xijki VjVi and K'^i = ^ Xij^i Vjvl . (A14) 

Computing the real matrices defined in Eq. (|A^), we arrive at the result 

A = Re{i2'^ + A + K') +ReK , (A15) 
B = Re{i2'^ + A + K') -ReK, (A16) 
C = -lm{ij,^ + A + K') -ImK . (A17) 

These matrices determine the 2nH x 2nH mass matrix Ai^ of the neutral scalars. Equa- 
tion (^) reads 

(/i^ + A + K') b + Kb* = mlb . (A18) 



3. The Goldstone bosons 

The Goldstone bosons corresponding to the longitudinal modes of the W and Z vector 
bosons are given, respectively, by 

aw = -{.Vi,V2,. . . .Vnaf and bz = - {Vi,V2, . . . .Vnaf , (A19) 

V V 

where 



/I |2 , I |2 , ,1 |2 'imw r\nn\ 

^=Vl^il +l^2| H ^IvuhI =—^- (A20) 
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They correspond to zero eigenvalues of the mass matrices Ai\ and Ail, respectively. Indeed, 
making the replacement 0fc — > (O|0fc|O) = {0,Vk/ V^)^ in the scalar potential of Eq. ( |A11| ) 
and enforcing the condition that this is a stability point of V, one obtains 



_d_ 

dv* 



V{{0 10,1 0)) = i + J2\^,HvXvi] =lj2i^- 

V j j,k,l / j 



)^^v, = 0. (A21) 



Thus, Ai^aw = 0. Furthermore, inserting bz into Eq. ( |A18| ) one obtains 



(/i^ + A + K') bz + Kb*z = i (K'aw - Ka^ 



w I 



(A22) 



where we have used bz = iaw, (/W^ + A) aw = 0, and the definitions of i^T and K' in Eq. (|A14D . 
Thus, rriayy = and rrih^ = O.Q 



4. Feynman rules for some gauge vertices 

The covariant derivative of the scalar doublets is (we use the notation of Ref. ||16|| ) 



D 



9 



W' 







9 
2r 



C,,, 5,,, 







+ie 









(A23) 



The weak mixing angle 9^, the positron charge e, and the SU (2) gauge coupling g are related 
hj e = gsw, where Cy^ = cosOy^ and = sin6'^. 

The covariant derivative in Eq. (|A23|) leads to W^S^S^ couplings given by 



£ = + , ^ ^ atfo {S- d'^S', - S', d'^S-) + H.c. 



(A24) 



a,b 



where the sums over a and b include the Goldstone bosons. The Feynman rules for the 
vertices W^^S'S^ and W~S^S^ are, therefore. 



i^a'^bipb-pa)'' and i^b'^a{pa - Pb) 



9 



(A25) 



respectively. In Eq. (|A25|) , ph and Pa are the incoming momenta of and of S^, respectively. 
/^From the covariant derivative ( |A23|) one also derives the coupling 

(A26) 



£ = •■■ + gmwW;^W-' ^ (a{^6) 5° . 

b^bz 



^ After addition of the gauge- fixing terms to the Lagrangian, we eventually have = mw- 

Since we use the unitary gauge for the Z boson, the neutral Goldstone boson S^^ does not appear 
in our calculation. 
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Note that a\yb is real for all b ^ bz, because of the second Eq. (|A6|) and bz = iaw (see 
Subsection |A3| ). 

The covariant derivative in Eq. (|A23| ) also yields the H^^S'j7 and W^S^Z couplings, 
given by 

C = mw f eA^ + ^ zA {W^G- + , (A27) 

where = S^^ are the charged Goldstone bosons. We emphasize that the charged Gold- 
stone bosons are the only charged scalars which have a vertex with W^'j and with W^Z. 
We also need the couplings of the photon and the Z boson to the charged scalars: 

C = --- + t (-eA- + g Z^^ i^ad.S^ - S^d.S^) . (A28) 

Here the sum includes the charged Goldstone boson. 

For the three-gauge-boson vertices and W^W^Z see, for instance, Ref. [Q. 



5. Feynman rules for some scalar vertices 

Let us consider the contributions to C = ■ ■ ■ — V with one positively-charged scalar, one 
negatively- charged scalar, and neutral scalars: 

= • • ■ - E E ^^^^^ <S-a'^S^, {vl + bis',) {vi + bX) ■ (A29) 

i,j,k,l a,a',b,b' 

We pick from this expression the vertex S^,S^: 

C = Caa'bS;,S^,S^ , with Caa'b = Qafe = E ^ioklOla'^ {vlk + Vlbl) . (A30) 

a,a',b 



In the particular case a' = aw, with aw given by Eq. ([A19| ), we obtain 



C=...--^J2Y1 (^^^' + ^^^aS', + H.c. (A31) 

zniw '—7, 

a,b i,j,K,l 

With the definitions of Eq. ( |A14D and using Eq. ( |A18D , we perform the simplification 



o:iv,\m Kk + hlvi) = Y < (K'b + Kb*), = {ml - ml) a^b . (A32) 

i,j,k,l i 

We thus arrive at the interaction 

C = --- + -^J2^^b{ml-ml) G^S-S', + H.c. , (A33) 



w 



a,b 



which is the counterpart of Eq. ( [A24D in gauges. Therefore, we have 
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Caawb = i^b - ^l) and Cay^ab = 7r~ i^b " ^l) ) (^34) 

and the Feynman rules for the vertices G~^S~S^ and G~S^S^ are 

i — - — a^b (ml — ml) and i — - — b^a (ml — ml) , (A35) 
2mw 2mw 

respectively. 

In Eq. ( |A33| ) we may consider the particular case a = aw- Remember that a\yb is real 
for all b ^ bz- Using = (see Subsection |A 3|) , we obtain 

Cayva^b = ml (a\^b^ , (A36) 

and the corresponding Feynman rule for the vertex G^G^ S^: 

^ ml(a\yb). (A37) 



2m 



w 



APPENDIX B: THE FERMION SELF-ENERGY AND ITS RENORMALIZATION 

1. On-shell renormalization conditions 

Most of the material in this Subsection can be found, e.g. in Refs. ||l^j2l|. We include 
it in order to make our paper self-contained, since the renormalization procedure outlined 
here is used in Section 

Let us consider a theory with n Dirac fermions, e.g. charged leptons. We use a matrix 
notation whenever possible. The unit matrix is not distinguished from the number 1. The 
one-loop fermion self-energy is 

- ^S(p) = -z y [Al {p^) IL + Ar + Bl [p") IL + Br . (Bl) 

The quantities Al^r and Bl^r are n x n matrices. They constitute the result of the one- 
loop calculation. We assume that there are no absorptive parts in the one-loop self-energy 
diagrams, then 

{Al);^ = {Al)^^, (Ah);^ = (A^)^^ , and (5^);^ = (5^)^^ . (B2) 

We shall also assume that the fermion masses are non-degenerate. 

We want to study the one-loop renormalization of the fermion fields. The bare chiral 
fields are denoted iaL and iaR, and the renormalized ones and i^^R-. 

h,R = v^ilR = (l + ^ ^^-«) ^Ir ■ (B3) 

The zl^r are n x n matrices. The diagonal matrix of the bare fermion masses m^ is Mi = 
diag(mi, . . . ,mn). We write = -|- 6M£, where is the diagonal matrix of the 
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renormalized masses, and SMi = diag {Smi, . . . , SrUn)- Then, the renormalized fermion self- 
energy is 



-i (^5M^ + i MeZL + ^ zl^Mi^ 7^ - i (^5Mi + ^ z[M^ + ^ Mez^^ 7^ . (B4) 

The terms in beyond those in originate in the counter-terms for the masses and 

for the wave functions. 

Up to now the renormahzed self-energy has no precise meaning. Now we fix the meaning 
of this notion by requiring that fulfills the following conditionsifl 



9 9 



Ur 



, W a, j3; 



Condition 1: 
Condition 2: lim [Sf{p)]^^ - m„) 

where the fermion propagator is given by 

Sf{p) 



Ur 



y a . 



1 



(B5) 
(B6) 

(B7) 



We are using the notation Ua = u{p, rua, s), where the four-spinor is to be taken with mass 
TUa, three-momentum p, and polarization s; it satisfies '^Ua = rriaUa- The above equations 
are the conditions for on-shell self-energy (or propagator) renormalization |1^,21,24-26|, and 
they include rotating back the renormalized fermion fields into the physical basis. Condition 
2 fixes the residuum of the propagator Sp{p) at the pole p"^ = to be 1. 

Exploiting first Condition 1, i.e. Eq. (|B5|), it is easy to see that it leads to the following 
relations, which originate in the coefficients of 7ij and 7^, respectively: 



1 1 

m„Ai (ml) + Br (ml) - -ma{zL)f3a + -^p{zR)/3a + SmaSpa = 0, 



1 1 



'/3q 



(B8) 
(B9) 



These equations may be solved to give |15|J21|1 , for a ^ P, 



1 . . _ ^l^L ("^a)/3a + "^a"^/3^ij ("^a),3a + ^P^L ("^a)/3a + "^"^^ ("^a)/3a 

2 [ZL)f3a - 



— 



1 , , ^l^R i^Dpa + mampAL (m2 + rripBR (m^ + m„Si (m^ 

2 {ZR)l3a - 



— 



(BIO) 
(Bll) 



and, for a = P, 



^Condition 1 is stated in Ref. 11151. Both conditions are stated in Ref. pTl. 
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-2 5m^ = m« [Al {ml) + An {ml) J + {ml) + Br {ml) , (B12) 
m„ [ZL - zn)^^ = m„ [A^ (m^)^^ - {ml) J - B^ (m^)^^ + Bn (m^)^^ . (B13) 

The right-hand side of Eq. (pT^) is real because of Eqs. Equations ([BTQ|)-([BT2D fix 

the mass counterterms and the off-diagonal wave-function counterterms. 

With Eqs. (|B2|) , the Zl of Eq. ( |B10| ) and the zji of Eq. (pil| ) agree with the corresponding 



self-energy renormalization for a P in Ref. [|1^], after one specializes the functions B^ ji 
to the forms used there. 

Only the difference {z^ — ZRjaa is determined by Eq. ( P13|) . In order to fix {zi)aa and 
{,ZR)aa separately, one has to invoke Condition 2. Equation ( p6D is equivalent to 



hm -^^S^Jp)w, = 0, Va, (B14) 

which is better suited for the further procedure. Firstly, we expand the matrix functions 
Al,r{j)^) and Bl,r{p^) in E(p) around p"^ = ml; for instance, 

A, {p') = A, {ml) + {p^ - ml) A', {ml) + ^ (p^ - m^) ^ A^ {ml) + ■ • • . (B15) 
Secondly, we use Eqs. (p^ ) and (|B^ ) to obtain 

Kaip) = [^L {ml) - RezL]^^ {ilL - m^7R) 

+ [^R {^D - Re zr] (^7^ - m„7L) 

+ {v" - O [A'l {ml)^iL + A'n {ml) 
{ml)i, + B'^{ml)^R]^^ 

+ \ {P' - <f [K {ml) ^7L + K) j/7« 

+Bl{ml)^, + B'l,{ml)^R]^^ + ..., (B16) 

where the dots indicate higher orders in p'^ — ml- It is obvious from Eq. ( [B16| ) that, when 
= ml, E^^(p)Ma = 0. This agrees with Condition 1 in Eq. (pSj). Condition 2, in 
Eq. ( |B14|) , is nonetheless non-trivial because of the presence of the denominator |( — m^, 



which also vanishes when applied to Ua- Thirdly, we use 



' ' L - majR) = (7l^ + ma'jR) -1—^ , (B17) 



f — m^ p + m 

together with the analogous relation for opposite chiralities, to obtain 

Kaip) = [^L {ml) - Re Zl] {-iii + "^a7i?) 



+ [Ar {ml) - Rez^]^^ (7i?^ + ^,7^) j^^^ 

+ {i> + ma) [A'l {ml) i)-iL + A^R (m^) i>-iR 
+B'l {ml) ^L + B'^ {OlR]aa 

+^^^ (i^ + ma) [A'i {ml) ^7L + A'A K) ^7k 

(m^) 7^ + B'l, {ml) 7i?] + • (B18) 
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Therefore 



[Al (ml) + An (ml) - Re {zl + zr)1 



+ml [A'^ {ml) + A'n (m^)]^^ + m„ [B'^ {ml) + B'^ K)]^^ 
+0(p2-m2)}M„. (B19) 



Finally performing the limit ml to apply Condition 2, i.e. Eq. ( |B14| ), we obtain |ST 

Re {zl + zr)^^ = [Al {ml) + Ar {ml) + 2ml [^l {<) + K)] 

+2m, [B'^ {ml) + B'^ {ml)] ] . (B20) 

Equations ( pi3|) and (P20| ) together fix the real parts of {zL)aa and of {zR)aa- Concerning 
the imaginary parts of those quantities, Eq. (P13|) fixes their difference, Ivci^zl — ZR)aa, 
whereas their sum remains undetermined. This fact reflects the freedom of redefining the 
fields iaL and iaR by transforming them with the same phase factor. 



2. The equivalence of two different procedures 



In Sections 0, |VI1| , and |VIII| we have considered, respectively, the vertices rfi'-f, rfiZ, 



and TfiS^. When computing those flavor-changing vertices we have not invoked any renor- 
malization procedure, but we have considered self-energy transitions r — /i in the external 
fermion legs. In the following we show that, instead of adding the self-energy transitions 
in the external fermion legs to the one-loop vertex, one may equivalently apply the on-shell 
renormalization prescription to the external fermion legs. We explicitly work out this equiv- 
alence, for an arbitrary fermion self-energy, in the case a scalar vertex (a Yukawa coupling). 
However, it will become clear that the nature of the vertex is irrelevant for this equiva- 
lence, and therefore our considerations are of general validity. In the case of the one-loop 
flavor- changing Z vertex in the Standard Model, this equivalence was explicitly derived in 
Ref. 



15 



Let us first see the effect of the on-shell wave-function renormalization on the Yukawa 
interaction of a real scalar S^, 



WL + hrliR 



(B21) 



where Tq is the n x n coupling matrix of unrenormalized coupling constants. Defining the 
renormalized scalar field via S'^ = \/ZsS^ = (1 + zs/'2)S'^, and denoting the renormahzed 
coupling matrix by F = Fq — 5r, the Lagrangian of Eq. ( P21|) , written in renormalized 
quantities, is 



Cy {S') 



+ 



.R^T + ST 
'Ft 



+ ^zsT^ + ^T^ZR + \ziT^ 



^R 



(B22) 
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The quantity 5T is determined by the renormahzation condition for T. 

We now consider the part of the vertex ~^ ^isS^ given by the wave-function renormal- 
ization constants -2l,_r. First we consider those terms which have the couphng matrices F 
and F"'' to the left of the factors z; only those terms can give a fermion self-energy on the 
leg. For any 7 7^ a, using the Lagrangian in Eq. (|B22| ) and taking into account Eqs. ( [B10|) 
and (|BTT| ) 



we obtain (see also Ref. |[T5[] ) 



F/3^ (zl)^^ IL + (F^) {zr)^^ -in 



Ur 



Up 



X 



h — 



(B23) 



where the four-spinors Ua and are to be taken at four-momenta pi and p2, respectively, 



with PI 



m: 



and P2 



Equation ( P23|) holds for every 7 7^ a and demonstrates that 



the z-terms taken into account in the left-hand side of Eq. ( [B23|) exactly reproduce the 
self-energy insertion in the leg ^a of the tree-level scalar vertex. 

Equation ( |B23| ) actually does not depend on the assumption that we are dealing with 
a scalar vertex. All operations needed in order to derive Eq. ( P23|) took place to the right 
of the matrices F and F^, and they involved only the spinor Ua- Thus, our scalar vertex 
might be replaced by any other vertex with different Dirac structure, and the result would 
be analogous. 

For the terms where F and F"!" are to the right of the wave-function renormahzation 



constants, we have to consider 



[zl,r)1i3, which means that the matrices Al r and 



Bl r, as functions of p^, must be evaluated at m^. This contrasts with Eq. ( [B23| ) where they 
had to be taken at m^; in Eq. (|B23| ) the fermion self-energy appears with pi. Furthermore, 
we now must use Eq. (|B2|) . Having this in mind, for any 7 7^ /? we arrive at 



Ur 



Up [-iS(p2)] 



Pi, 



F^„7L + (F^) 



70 



Ur 



(B24) 



In conclusion, in the example of a simple scalar vertex and for a ^ ^ (3 one can see 
the equivalence [O] of the following two procedures for calculating the transition amplitude 



1. On-shell renormahzation of the propagator in order to calculate the wave-function 
renormahzation constants Zl r, and then taking them into account at the vertex; 



2. Adding the two self-energy contributions to the vertex, without any renormahzation. 
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